Zero-bias anomalies of point contact resistance due to adiabatic electron 
renormalization of dynamical defects 
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We study effect of the adiabatic electron renormalization on the parameters of the dynamical 
defects in the ballistic metallic point contact. The upper energy states of the "dressed" defect are 
shown to give a smaller contribution to a resistance of the contact than the lower energy ones. This 
holds both for the "classical" renormalization related to defect coupling with average local electron 
density and for the "mesoscopic" renormalization caused by the mesoscopic fluctuations of electronic 
density the dynamical defects are coupled with. In the case of mesoscopic renormalization one may 
treat the dynamical defect as coupled with Friedel oscillations originated by the other defects, both 
C^") ' static and mobile. Such coupling lifts the energy degeneracy of the states of the dynamical defects 

giving different mesoscopic contribution to resistance, and provides a new model for the fluctuator 
as for the object originated by the electronic mesoscopic disorder rather than by the structural one. 
D ' The correlation between the defect energy and the defect contribution to the resistance leads to 

5^ , zero-temperature and zero-bias anomalies of the point contact resistance. 

■ A comparison of these anomalies with those predicted by the Two Channel Kondo Model (TCKM) 

is made. It is shown, that although the proposed model is based on a completely different from 
TCKM physical background, it leads to a zero-bias anomalies of the point contact resistance, which 
are qualitatively similar to TCKM predictions. 
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Recent advances in nanofabrication technology have 
made it possible to visualize single defects with inter- 
nal degrees of freedom - "fluctuators" which lead 
to a "telegraph" resistance noise of nanometer scale sys- 
tems. In metals these defects are believed || |?J to be 
some structural defects, which at low temperatures are 
OS • seen as the well-known two-level tunneling states (TLS) 
d|. TLSs are the objects typical for strongly disordered 
amorphous solids || that switch by tunneling between 
■ their two possible configurations. Although the micro- 
scopic nature of the fluctuators stays unclear (especially 
for ballistic devices made of pure metals), the experi- 
ments allow to study various phenomenological parame- 
Q , ters of these objects, in particular, the interlevel spacing. 

Kondo H was the first, who has pointed out that 
• i-H , in metals the parameters of the dynamical defects are 
strongly renormalized by electrons. One can discrimi- 
Jh I nate between adiabatic electron "dressing" , which is re- 
lated to a static electron response on the defect poten- 
tial, and non-adiabatic one, which affects the tunneling 
process and leads to a renormalization of TLS tunneling 
matrix clement ( "dissipative tunneling"). It is this non- 
adiabatic effect that has attracted most attention (see 
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e.g. Ref. [Q) due to its evident importance for the de- 
fect dynamics. As for the adiabatic renormalization of 
the defect parameters, it is customary to include it into 
the bare values. This procedure is usually justified by 
the fact that the adiabatic effects are related to a re- 
sponce of the whole electron systems while only a small 
strip of electron energies, close to the Fermi surface, is re- 
sponsible for transport properties and sensitive to exter- 
nal factors like temperature or fields applied. However, 
as it was first demonstrated in Refs. the adiabatic 

"dressing" of the fluctuator, in particular, the adiabatic 
renormalization of the fluctuator energy splitting E, can 
be important and depends on the state of the electron 
system (e.g. on superconducting properties). Further- 
more, very recently the surprizing " magnetic tuning" of 
the TLS interlevel spacing observed 11 for TLSs in Bi 



nanoconstrictions was explained 12 as a direct result of 
the adiabatic renormalization of TLS parameters by elec- 
trons, which states are affected strongly by the magnetic 
field. 

The purpose of the present paper is to show that the 
adiabatically renormalized energy of the fluctuator state 
correlates with the fluctuator contribution to the resis- 
tance. Namely, the conductance is larger for the fluctu- 
ator on its upper level. This fact causes, in particular, 
Kondo-like zero-temperature and zero-bias anomalies in 
differential conductance of the metallic point contacts. 
Indeed, an increase of temperature leads, in average, to 
an increase of occupation numbers of the fluctuator upper 
levels. The abovementioned correlation causes a corre- 
sponding conductance increase with temperature, which 
imitates the Kondo-like behavior. The same holds for an 
applied bias increase. 



II. ADIABATIC RENORMALIZATION OF THE 
PARAMETERS OF THE DYNAMICAL DEFECTS 
IN METALS 

Two mechanisms of the adiabatic renormalization can 
be considered. First one is due to possible difference of 
electron-fluctuator coupling potentials V for the two of 
the fluctuator states ^ |^ (2) |) and was studied 

in Ref. pjji], L2[ . We will show that for this mechanism 
(which will be referred to as "classical" one) the above- 
mentioned correlation is due to the fact that expressions 
for the conductance and for the electron contribution to 
defect energy include the same strength of the electron- 
defect coupling. 

Second mechanism of adiabatic renormalization of dy- 
namical defects parameters was suggested by Altshulcr 
and Spivak |13| . It implies mesoscopic electron den- 
sity fluctuations, which lead to a difference, even for 
\V^ \ = \V^ \ , of electron-fluctuator coupling strengths 
for different fluctuator states due to their spatial reso- 
lution. The same correlation does occur for this "meso- 
scopic" contribution, as for the "classical" one. We will 
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show it for the experimentally important situation of bal- 
listic point contact, where a description in terms of "lo- 
cal" interference, which involves finite number of scatter- 
ers, is possible. 



A. Renormalization due to difference in 
defect-electron coupling potentials 

Let us start from the "classical" effect. TLS energy 
splitting E depends ]tJ on the TLS tunneling parameter 
Aq and on the TLS asymmetry A: 



E=^/At + A2. (1) 

Interaction of TLS with conduction electrons leads to the 
renormalization of both TLS asymmetry ^|||[l2| and of 
tunneling parameter || . For the asymmetrical TLSs with 
the large enough barrier, A 3> Ao, the renormalization 
of TLS asymmetry gives the major contribution to E. As 
a result 12 1, the electron-TLS coupling leads to the 



renormalization of TLS "bare" energy splitting E which 
has a form E — > E + E e i, where E e i — E e i t 2 — £eZ.i, 

= £ I y;"'!' gaWM .-fifl!. (2) 

Here Vq is the electron-TLS coupling constant for i- 
th TLS configuration. The electron bandwidth is as- 
sumed to be of order of £p, and the electron distribution 
/(e) to have an equilibrium Fermi form /(e) = /o(e) = 
{exp[(e — Ef)/T] + l} -1 . Applicability of the second- 
order perturbation theory approximation is justified for 
Eq. (||) when V^f* <£f. In contrast to Kondo-like cor- 
rections, Eq. (0) depends weakly on temperature. The 
renormalization E e i is due to the difference in the values 
of the total electron energy, renormalized by a presence 
of TLS, for the TLS in states 1 and 2, respectively. It is 
important, that E e i is not, due to adiabaticity, sensitive 
to the details of inter-state transition mechanism. Thus 
the problem is reduced to estimates of energies E e i^ cor- 
responding to different configurations of the defect which 
for this case can be considered as its independent realiza- 
tions. Therefore, Eq. (||) holds both at low temperatures, 
when the fluctuator transitions are due to tunneling, and 
at higher temperatures when thermal activation domi- 
nates [flQ- 

According to Eq.(||) the energy E is lowered from the 
"bare" value. The lowering is larger the stronger is the 
coupling of the state with the electrons. On the other 
hand, a presence of the defect in the state i inside a 
ballistic point contact with a characteristic size d causes 
a reduction of the contact conductance Q 

1 \ ( \f 2 
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where Xf is the Fermi wavelength, and G ~ 
(e 2 /h)(d 2 JX F ) is the Sharvin conductance. Making use of 
the Eq. (g) one obtains 5Gi w (e 2 /h)(E e l,i/ £f)- There- 
fore, if the defect asymmetry is completely controlled by 
the conduction electrons, the conductance, which corre- 
sponds to the defect in its lower state (larger absolute 
value of the electron contribution to the defect energy) 
is smaller than that for the defect in the upper state. 

B. Interference contribution to the renormalization 

Let us consider now the electronic interference. For 
adiabatic effects different states of "active" defects, fiuc- 
tuators, can be considered as independent realizations. 
Thus we can choose some configuration of some active 
defect as a "reference" scatterer I and consider its prop- 
erties in the presence of "background" scatterers. We 
will analyze both the interference contribution to con- 
ductance due to the defect i and the "mesoscopic" renor- 
malization of the energy of this defect, which is equal to 
a change of the electron system energy due to a presence 
of interference pattern involving defect i. 

The ballistic point contact contains finite number of 
scatterers. Therefore the interference contribution to the 
contact conductance is provided by a local interference 
(involving trajectories with small number of scatterers) 
rather than by a global one (which leads to well-known 
UCF |u|). This "local" interference contribution has 
been to some extent analyzed in Ref. |17| . In what fol- 
lows for simplicity we will restrict ourselves mainly to 
the interference patterns involving only pairs of scatter- 
ers. However, as it is shown in Appendix la, our results 
can be generalized for the case of arbitrary number of 
scatterers. 

As it is shown in Appendix la [Eq.(j27])], the contri- 
bution to the conductance due to a pair of scatterers, 
namely "reference" scatterer i and the "background" m, 
is: 

^e 2 \V \ 2 ( \ F V 

Here £j m = cos(2fci?i?i m ) and Ri m = Ri — R m is the 
vector, which connects two scatterers. For simplicity we 
assumed that the scattering potentials for all scatterers, 
both "active" and "passive" ones, depend only on the co- 
ordinate R; of the scatterer: Vikk' = Vq exp[i(k — k')RJ. 

Let us now find a contribution to the energy of "refer- 
ence" scatterer i, which represents one state of some fluc- 
tuator, due to its being involved in electron interference 
along with another scatterer m. Following the scheme 
implied by Eq. (|^), it is given by a renormalization of 
electron system energy due to this pair of scatterers. In 
second-order perturbation theory approximation we ob- 
tain, collecting all terms proportional to V^kk' Knkk' • 
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E el>im = ReV|l/ | 2 exp[*(k-k')R™] /(gk) /(£k ' } ■ 

(5) 

A straightforward calculation for a spherical Fermi sur- 
face and zero temperatures gives [see Appendix la, Eq. 

TP - A C A ~ I V I 2 A F / fi \ 

The obtained renormalization is due to interaction of de- 
fect m with the Friedel oscillation of electron density orig- 
inated by the defect i. 

Both 8G im , Eq. and E e i >im> Eq. (|), arc pro- 

portional to the same phase factor £j m = cos(2kp Ri m ). 
Correspondingly, 

[6G m /(e 2 /h)} ^Xf_ >Q (?) 

(SE e l y i m /eF) Rim 

and thus, in analogy with the "classical" effect, the larger 
is the energy of a configuration, the larger is the contact 
conductance. As it can be shown (see Appendix la) the 
proportionality to the same phase factor and, therefore, 
Eq.0, holds not only for pairs of scatterers but for arbi- 
trary number of scatterers as well. 

To estimate a total interference contribution due to 
the defect i both to the conductance, SGi, and to the 
defect energy, E e i^, one must sum over all "background" 
scatterers m. Due to a random distribution of Ri m this 
results in some mesoscopic fluctuations for both quanti- 
ties with respect to realizations of the system. However 
the fact that both quantities are linearly related to the 
same set of random factors & m leads to the correlation 
between them, namely 

(5G i E elti } = CSGE7i, (8) 



where x = J (x 2 ), and ( } denotes the ensemble 

average. More detailed argumentation of this is given in 
Appendix 2. For a given value of E e ii = E one has 

(SGi) E = C6G^^. (9) 

Keeping in mind factors \/R 2 m (for the conductance) 
and l/i?f m (for the energy) one may suggest the main 
contribution to both these quantities to stand from the 
nearest neighbors. In this case both energy renormaliza- 
tion and contribution to resistance would be related to 
a few neigbouring defects and thus could be estimated 
by EqspL with Ri m of the order of most probable in- 

— 1/3 

terdefect distance N t ' (where Ni is the background 
defect concentration). Note that if we took an average 
over all possible realizations of the background scatterers 
we would met the problem with singular behavior of the 
averaged quantities when R im — > 0. It means that the 
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average is mainly controlled by (rare) realizations corre- 
sponding to very close neiubouring background defect, 
and one has Ri m > kp 1 , that would give \SG\ ~ e 2 /h, 
E e i t i ~ Ep I (kpl) 1 / 2 . However, in the case of small ballis- 
tic contact we deal with some given contact realization, 
so that one deals with the most probable rather than 
with average quantity (compare with |fl5||). 

For the ballistic point contact the number of defects in 
the contact is small, so that the dominant contribution 
is expected to be from the trajectories, which involve a 
boundary of the contact. Assuming the contact to be a 
short channel with a length ss d, this boundary may be 
considered as an array of scatterers at a distances m d 
from the defect i with a total number « (d/Xp) 2 . For 
this case 

w -(t)p (¥)'■ <"» 

Taking values of d typical for nanofabricated ballistic 
point contacts, d w 5 — lOnm one gets \SG\ R* (0.05 — 
0.1)e 2 /h, E ehi w 30 - lOO-K". 

The mesoscopic interference renormalization has some 
special features as compared with the "classical" one. 
First, mesoscopic disorder lifts energy degeneracy of the 
defect states, which have different spatial positions. Thus 
it causes a formation of fluctuators from otherwise sym- 
metric defect configurations (to say, interstitials, which 
have symmetrical lattice positions). 

Then, in this case one expects the temperature and 
bias behavior of the resistance to depend on an external 
magnetic field, which affects the electron interference (see 
e.g. Ref. @). 

In addition, the interference contribution both to the 
defect energy and to the conductance depends on the 
electron distribution. The finite applied bias makes this 
strongly nonequilibrium, which at high enough biases 
causes a "direct" effect of bias on both quantities. 

The physical picture of the mesoscopic renormalization 
is much richer than that provided by "classical" one, and 
it is this mechanism that we will concentrate on in the 
rest of the paper. 

III. ZERO-BIAS AND ZERO-TEMPERATURE 
RESISTANCE ANOMALIES 

Let us consider defect i which occupies either of the 
two neighboring positions, 1 and 2, with close energies. 
For simplicity we will assume that the energy asymmetry 
of these defect states is completely determined by the 
electron renormalization. This object is a sort of two- 
level fiuctuator originated from the electronic disorder 
rather than from the lattice one. It is important, that due 
to the correlation discussed above, the upper state of such 
fiuctuator, which corresponds to a defect position with 
higher energy, gives a smaller contribution to the contact 
resistance. A conductance increase, which accompanies 
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a transition from the lower to the upper level of such 
fluctuator, 8G(v, = SGu^x — 8Gu\q, is, according to Eq. 
(^), scaled with the energy asymmetry Eu) = -E e z.(i)i — 
E e i,(i)2- Here index i now denotes the fluctuator. 

A summation over all fluctuators gives their total con- 
tribution to the average contact conductance: 

AG = J2SG (i) n (i) = J (5G {l) ) E P{E)n{E)dE. (11) 

i 

Here — n(Eu\) is the i-th fluctuator upper level oc- 
cupation number and P(E) a density of states given by 
statistical properties of £ e i,i- For the mesoscopic sys- 
tem it is reasonable to take the values of E e i^ for the 
neighboring defect positions as statistically independent. 
In this case P{E) is approximately constant at small 
E <C E e i. Making use of Eq.(^|) and of the expression 
n(E) = [1 + exp(E/T)]~ 1 one has at small temperatures 
T <C E e i the conductance enhancement 

AG oc T , (12) 

where (3 — 2. 

For some defects, like light interstitials or some defect 
complexes, the probabilities for defect hopping between 
spatially symmetric positions are relatively high jl8| . For 
these "delocalized" defects the effect of electronic dis- 
order provides a many-site "potential relief" instead of 
two-site fluctuator picture. Assuming that any site can be 
occupied by only one defect one deals with "Fermi- type" 
statistics; so at T — > the sites with lowest energies are 
occupied by the mobile defects while those with energies 
higher than the "Fermi level" are free. At finite tem- 
perature Eq.|ll] can be applied where the site occupation 
number has again the form n(E) = [1 + cxp(E /T)]^ 1 if 
one takes the "Fermi level" as the origin for the energy 
E. In this case the total number of available sites is much 
larger than the number of defects TV, and for finite tem- 
peratures the " Boltzmann-type" statistics holds rather 
than the " Fermi- type" : 

N eM-E/T) 
n{E) = JP(E)cM-E/T)dE (13) 

In this case a change of T does not affect the number of 
rearranged defects (because any of them can change its 
energy) and leads only to a change of the average defect 
energy. As a result, in this case we have in Eq.(]12|) /? = 1 
independently of the form of the density of states P(E). 

Let us turn now to the effect of finite bias eV ^ T . 
For TLS it was first considered in Refs. |5|JlS|]. It was 
shown that for low-energy TLS with small enough en- 
ergy splitting E, for which the coupling with electrons 
dominates EG], the TLS occupation numbers are sensi- 
tive to the electron distribution. For the contact region 
this is strongly nonequilibrium and for a central point of 
symmetric contact has a form 

/(k) = 9{k x )f (e k + eV/2) + 9(-k x )f (s k - eV/2), 

(14) 
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where OX is the main axis of the contact and 8(x) the 
theta function. The "energy width" of this distribution, 
eV, plays a role of the effective temperature. In partic- 
ular, the upper levels of the TLS are empty if eV < E, 
while for eV ^> E the occupation numbers of TLS levels 
are almost equal HQ and n(E) « (1/2) [1 - (E/eV)]. 

For larger E the coupling with phonons becomes 
important @,[l4| due to rapid increase of density of 
states for the actual phonons with E increase. For 
the two-state case the fluctuator relaxation rate due 
to electron-assisted tunneling is W e i (E, V) ~ 

(\V -\/eF) 2 [(eV - E)/h]T, where V ~ = (V (2) - V {1) )/2. 
For the phonon-assisted process (?]] W p h(E,T) k 
(A 2 E 2 /£e 3 D )(E/H)cth(E/2T)T. Here A is the fluctu- 
ator dcformational potential, £ and Qd are the atomic 
and Debye energies, respectively. T = exp(— A), where A 
is the tunneling constant. 

Let us define the characteristic energy E* for which 
W el {E*,(eV - E*) « E*) = W ph (E*,T = 0); for the 
reasonable values of parameters (see e.g. Ref. |p0f ) E* 
is expected to be k 1 — 3A". For eV 3> E* a prob- 
ability of electron-assisted tunneling to the upper level 
W e i exceeds the probability of phonon-assisted decay 
of the upper level W p h up to some threshold energy 
E = E t h = E*{eV/ E*) 1 / 3 at which an increase of the 
electron phase volume with bias (oc eV) is compensated 
by a corresponding increase of phonon phase volume 
(oc E 3 ). For the crude estimates let us take the occu- 
pation numbers n{E) oc 6{E t h — E). Now making use of 
Eq. ([ll]) and assuming density of states P{E) constant, 
one obtains the following interpolation formula for the 
interference contribution to the conductance: 



AG 



eV\ 1/3 
E | j r \ T 



(15) 



with /3 = 2. 

The same considerations can be applied for the case 
of "delocalized" defects. Although the probabilities W e \ 
and Wph for this case can differ from ones for the two- 
level fluctuators, the scaling W e i/W p h ~ eV/E 3 (relation 
between relevant electron and phonon phase volumes) 
holds for (eV > E > T) and thus Eq.[l| is valid, but 
with 13 = 1. 



IV. DIRECT EFFECT OF THE APPLIED BIAS 
ON THE FLUCTUATOR PARAMETERS 

When a large enough bias V is applied to the point 
contact one should take into account the non-equilibrium 
electron distribution in course of estimates Eqs.(||), (|^), 
(^). For the distribution given by Eq.(|l4|) one obtains 
(see Appendix lb) for AE e i^ m the phase factor 

AEei^m oc cos[2k F Ri m + tf>(V, R im , fejp)] cos(2AkR im ), 

(16) 
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and for 5G.; 



SG im oc cos[2k F R irn + 6(V, R im , k F )] cos(2AkR im ) 

86 1 



dV 2R{dAk/dV) 

where Afc = kpeV/ep. This is the "direct" effect of 
bias on the fiuctuators parameters in addition to tuning 
of fhictuator level occupation numbers. 

As seen, the first term in EqJT?] is correlated with the 
phase factor of Eq.[il| while the second is not and, thus, 
will sum out. As for the correlated cosine terms, the 
effect of bias initially ( at AkR <C 1) leads to their de- 
crease due to a decrease of the corresponding cosine fac- 
tors, while for AkR 3> 1 (when the factors are random 
with respect to parameter R im ) they are suppressed due 
to additional (with respect to the case V = 0) averaging 
over R im : 



sin[2k F R lm + 6(V, R lm ,k F )} sm(2AkR tm ), (17) 



(AE eh AG) oc (AkR) 



-1/2 



Actually we deal here with the well-known energy aver- 
aging effect suppressing any mesoscopic phenomena. 

It is important to note that these effects can lead to 
the resistance anomalies even if the defect structure is 
not rearranged in course of the bias application; the only 
condition is that the defects occupy the positions with 
the lowest energies available and thus with largest meso- 
scopic contribution to resistance SGi. For these defect 
configurations one gets in average SG < 0. Total inter- 
ference contribution to the conductance due to N defects 
is: 

AG - NF(V)6G, (18) 

where T ~ cos(AkR) for AkR < 1, and T ~ (AkR)- 1 ' 2 
for AkR 3> 1. The result we obtained is that the bias 
increase leads to a systematic conductance increase. It is 
interesting that in combination with the effects discussed 
for relatively small V — occupation of states with higher 
energies, this "direct" effect can form the configurations 
with smaller resistances than available for a simple tem- 
perature increase. Indeed, it can suppress (negative in 
average) mesoscopic contribution to conductance due to 
configurations with large enough energy gap between the 
available realizations which can not be rearranged at rela- 
tively small temperatures when the phonon contribution 
to resistance (obviously masking the effects in question) is 
still small. Note that as we have seen above the bias val- 
ues allowing the same occupation states of the defects and 
at the same time the same efficiency of electron-phonon 
processes as in equilibrium state with a temperature T is 
scaled with T as eV « T(T/E*) 3 . Thus the energy aver- 
aging effects can become to be pronounced for large-gap 
configurations when the filling of the upper level is still 
negligible. 

Certainly, the temperature increase can also lead to the 
energy averaging, but the necessary temperatures are too 
large and correspond to a significant phonon contribution 
to resistance. 
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V. DISCUSSION 



In this section we would like to make several remarks 
about the limitations and possible complications of our 
model. 

First, until now we considered the defect energy den- 
sity of states as constant. The limitation of the defect 
energy band leads to a saturation of AG(T) and AG(V) 
dependencies at T > T sat and V > V sa t, respectively. 
These quantities scale as eV sa t ~ T sat (T sat / E*) 3 , and 
AG sa t ~ NfSG, where Nf is a total number of fluctua- 
tors. As for the estimate for T sat , taking sf ~ 10 4 — 10 5 K 
and [d/Xp) ~ 50 and making use of Eq. ( |l0| ) we obtain 
T sat ~ E e i t i « 4 — 40K. Note that the "saturation" value 
AG sat corresponds to a random realization of different 
interference patterns involving the fluctuators, while the 
values of AG at lower temperatures correspond to prefer- 
able occupation of larger resistance states and thus are 
systematically smaller than typical for mesoscopic disor- 
der. 

Second, it is important that the picture discussed is 
sensitive to the external magnetic field. In particular, it 
is known that in homogeneous diffusive conductors the 
interference particle-particle channel is suppressed in the 
strong enough magnetic field H > &q/L^, where $o is 
the quantum of the magnetic flux and L c is a coherence 
length, instead of which in our case we should use a con- 
tact size d. This suppression reduces a magnitude of 
mesoscopic fluctuations nearly twice [jl3Lf[6|. Point con- 
tact is a strongly inhomogeneous system and the main 
contribution to the mesoscopic fluctuations is due to lo- 
cal interference. However, despite the fact, that the effect 
of magnetic field implies a contribution of configurations 
which involve more than two scatterers, this contribu- 
tion is relatively large due to rather high probability of 
the boundary scattering and lead to decrease of AG with 
field increase. 

Another important feature is related to a coupling be- 
tween different "active" defects, i and j, due to a depen- 
dence of the defect i energy on a position of defect j . This 
dependence is given by Eq.(|). For large enough concen- 
tration of "active" defects one may expect a formation of 
self-organized aggregates in the defect system (of spin- 
glass type). Indeed, the defect positions corresponding 
to maxima of the Friedel oscillations originated by the 
other defects become energetically preferrable, which in- 
troduces some "ordering" into the defect system. Thus 
a formation of " coherently scattering" aggregates can be 
possible leading to a significant enhancement of the inter- 
ference contribution to resistance. The increase of bias is 
expected first to suppress this contribution to the resis- 
tance in a way similar to discussed above. On the other 
hand, at higher biases the direct bias-induced decrease 
of "coupling potentials" Ea,ij (see Eq. (jig)) can de- 
stroy such aggregates, which can lead to sharp resistance 
changes. 

Finally, it is instructive to compare results of our model 
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with the two-channel Kondo model (TCKM) HQ, 
which also predicts zero-bias resistance anomalies of a 
non-magnetic nature. Despite these two models are 
based on completely different physical background, they 
predict qualitatively similar resistance behavior at low T 
and V: negative temperature and bias resistance coeffi- 
cients affected by a magnetic field. However, quantitative 
predictions of the two models differ. Our model does not 
predict a singular T-behavior at T — > - in contrast to 
TCKM. As for bias dependence, our model predicts sin- 
gular behavior, V 2 ^ 3 for biases V > 1 — 3mV (see Eq.(|l5|)) 
and saturation at smaller biases. This saturation can 
imitate the "restoration of the Fermi liquid behavior" 
predicted by TCKM. On the other hand, TCKM, being 
related to non-adiabatic effect, is relevant to the fluctu- 
ators of a rather special type (with a small asymmetry 
and large tunneling probability), while our predictions 
hold for any sort of mobile defect. Our model also pre- 
dicts special features at higher temperatures and biases; 
in particular, the saturation of the zero-bias anomaly at 
large V and T, and a principal possibility to reach larger 
values of conductance in course of bias increase with re- 
spect to ones available for temperature increase. 

VI. CONCLUSIONS 

To conclude, we predict a new mechanism of zero-bias 
resistance anomalies in metallic point contacts based on 
a found correlation between energies of defects with in- 
ternal degree of freedom and their contributions to re- 
sistance. The correlation lifts to some extent the "ran- 
dom" character of mesoscopic disorder and breaks the 
symmetry of the defect states with respect to the signs 
of the mesoscopic contribution to resistance. We suggest 
a model of a fluctuator related to a purely electronic dis- 
order, which provides a new insight into the nature of 
fluctuators in the perfect point contacts. 
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A. Calculation of the mesoscopic contribution to the 
conductance and to the defect energy at small 
applied biases. 



To calculate mesoscopic contributions both to conduc- 
tance and to energy of electron system due to the pres- 
ence of some finite number of scatterers we will make 
use of the "wave-optics" approach [M. The approach is 
based on perturbation theory in real space. Let us con- 
sider an electron, the wave function of which initially is 
a plane wave with the wave vector k. After n successive 
scattering events involving scatterers 1, n the electron 
wave function becomes: 

jn 

|1, ...,n >= V>l,...,n( r ) 



r-R n |-...-|R 2 -R 1 | 
x e lkRl exp {ik[\r - R„| + ... + |R 2 - Ri|]) (19) 

where Ri, ...R n are the positions of the scatterers. For 
the short-range scatterers the scattering amplitude / in 
the Born approximation takes the form: 



f = -^J drVM < 20) 

and is assumed the same for all scatterers. Contribution 
of the scatterers to the conductance is determined by 
their backscattering efficiency The interference contri- 
bution to the backscattering current due to trajectories 
involving scatterers 1, n and 1', n' is 

Si = 2^ (< V ' -' n ' |V|1 ' > +c - c ' ) - (21) 

To obtain a contribution to the conductance one should 
integrate this equation over r within some reference plane 
remote from the scatterers system. It is important that 
only position of last scatterers n and n' are relevant for 
this integration and one deals with a factor 

J d 2 pexp [ifc(| r -R„|-|r-R n ,|)]. 

Here p is a projection of r on the plane in question. Tak- 
ing the plane to be normal to R„ — R n / and expanding 
the exponent as [...] ~ |R„ — R n < | + p 2 \R n — R n < |/r 2 one 
gets the result of the integration in the form: 

exp(ifc|R„ - R„/|). 



ik\H n — R 



The next step is the integration over k directions. In its 
turn, this integration is relevant only to co-ordinates of 
the "first" scatterers in the chains that is to Ri and Hy, 
which enters the exponential factor exp[k(Ri — Ri')]- 
Correspondingly, the integration over cos 9 where 9 = 
Z(k, Ri — Ri') gives the factor 

i k\R 1 1 -R 1 ,\ eMlklRl - Rl ' l) - 
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Finally one arrives at the following estimate for the meso- 
scopic contribution to the conductance: 



5G /"+"' 
~G ~ k 2 a?R n + n ' 

x cos [<p(n, n', n' — 1, 1', 1; n', n,n — l, 1)] (22) 

where the phase <p is 

<p(n, n ,n — 1, 1 , 1; n, n — 1, 1) 

= /c(|Ri' — Ri| + |R„/ — R„| + |R„' — R„'_i| + ... + |R2' — Ri'| — |R n — R n -i| — ••■ — IR2 — Ri|)- (23) 

Here R is a typical interscatterer distance within the chains while d is the contact size appearing as a result of 
normalization of the backscattering efficiency on the incident electron flow. 

The picture discussed can be interpreted as a contribution to scattering due to presence of the scatterer n affecting 
the superposition of states formed by successive scattering by chains 1, ...,n — 1 and 1', ...,n'. The phase (p which 
after the integration over k directions is the phase difference for the paths n,n',n' — 1, 1', 1 and n,n — 1, 1, 
correspondingly. One should also note that in course of derivation of Eq.^2] we have taken into account that only 
those electrons with energies close to the Fermi energy contribute to the conductance, and used k — Uf in Eq. (p2|) . 

Now let us estimate the mesoscopic contribution to the electron energy due to the presence of the same system of 
the scatterers finally affecting the electron state in the position of scatterer n. In the lowest approximation one has 

SEei =< l',...,n'\V(r--R n )\l,...,n-l > (24) 

where V is a scattering potential assumed to be short-range: V = VoS(r). As a result of averaging over the direction 
of k we obtain 

SE el « -V Q J 



f, R n+n'+l 

x sin[k(p(n, n' , n' — 1, l', 1; n, n — 1, ...1)] (25) 

It is important that due to the same structure of the expression for SE e i as of one for the 5G the phase difference for 
the interference pattern ip is exactly the same. 

In order to obtain the interference correction to the energy of the whole electron system we should sum Eq.^5| over 
all occupied electronic states. For T = one has 

5E el 1 5~-rri / kdksm[k(p(n,n',n' -1,..., l',l;n,n -!,...!)] = 



n 2 jfn+n'+l 



Vo f n+n '- 1 k F , , 

•^2 gt +n'+2( n + n /) cos[<p(n,n,n - 1, 1 , l;n,n - 1, ...1)] 



(26) 



Thus we conclude that the mesoscopic contributions 
of the same system of scatterers to the conductance and 
to electron energy renormalization depend on the same 
phase factor and, therefore, are correlated. 

For the simplest case of 2 scatterers (positioned in Ri 
and Ri') site Ri plays at the same time a role of the site 
R„. The phase factor in this case is 

cos(fc|Ri' - Ri| + fc|Ri' - Ri|) = cos(2fc|Ri/ - Ri|) 

and, correspondingly, 

s -§ ~ ojkjw co < 2kR ^ SE * - -^jlr co < 2kR ) 

(27) 

which gives Eqs.(0) and (0). 
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B. Direct effect of bias on the interference 
contributions 



In this subsection we will study the "direct" effect of 
bias on the renormalization of the fluctuator energy and 
on the interference contribution to the conductance for 
the simplest case of 2 scatterers. 

We start with deriving E e i. For the step-like electron 
distribution given by E q.[14l one obtains for the integral 
over k (instead of Eq. fl26|)): 

rk F +Ak r k F -Ak 



^ f>K F -tL±K />K F —L±K 

- dk-sm{2kR + 6{V,R 7 k F )) + dk ■ sin(2kR + 6(V, R, k F )) = 

2 Jk F -Ak Jo 
rk F +Ak -. rk F +Ak 

/ dfc • ksm(2kR + 6{V,R,k F )) - - / dfc • k sin(2fci? + cj>{V, R, k F )) (28) 

JO 2 Jfcp-Afc 



Here A A; = eV/Tiv F . An additional ^-dependent phase 6 is related to a dependence of k on coordinate due to a 
presence of electric field ((hk) 2 /2m + ip(r) = e = const). Calculation of the integral gives the phase factor 

cos[2(fc F + Ak)R + 6(V, R, k F )] + sin(2fc F i? + 6) sin(2Afci?) = cos[2fc F i? + 6{V, R, k F )] cos(2Afci?) (29) 

instead of the factor cos(2fc F i?)obtained for V — * 0. 

In the same way we estimate the contribution to conductance at T = 0: 

r k F +Ak 



pK F -J-iiK 

Sloe / cos(2kR + 6{V,R,k F )), 

Jk F -Ak 



dl dAk 

8G= — oc -— [cos(2(fc F + Ak)R + 6(V, R, k F )) + cos[(2(/c F - Ak)R] + 6(V, R, k F ))} 

k F +Ak 

sin(2kR + 6(V,R 7 k F )) 



k F -Ak 



dV 

dAk 86 1 

- — — 2cos(2k F R + 6(V,R,k F ))cos(2AkR) + t^ — 2sm(2k F R + 6(V, R, k F )) sin(2Afci?) (30) 
oV oV 2R 

Taking into account that 6 ~ AkR ■ R/a, one sees that 
the second and the first terms in r.h.s. of Eq. (|30|) are of 
the same order provided that R/a w 1. However the first 
term completely correlates with the corresponding phase 



factor for the energy renormalization, Eq. (29), while the 
second term does not. 



VIII. APPENDIX 2 

For each "active" defect i the interference contribu- 
tion to the conductance, Gi, as well as to the energy, 
E e i t i contains summation over "background" scatterers 
m. Contribution of each scatterer m gives some phase 
factor £(Rj m ) = £im = £,mi which depends on the dis- 
tance of the scatterer m from the defect. Hence, one can 
rewrite the expressions for 8Gi and E e i^ in a form: 

SGi = ^""^ Gjm^im = (Gi,£i), (31) 
m 

E e i,i = ^^E im ^ im = (Ei,^j) (32) 
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Here we have introduced some " vector space" , where vec- 
tor £j contains the set of the corresponding phase factors, 
and vectors Ei and Gj contain the sets of the prefactors 
(given by Eqs. (3) and (5), correspondingly). For the en- 
semble of defects i the vector should be considered as 
random, while all components of Ei and Gi are positive. 
We may rewrite the vectors Ei, Gi as 

E - G - 

Ei = -gr-Eij Gi = -^-Gi (33) 

i^i tri 

where we have introduced the " norms" of the vectors Ei 
and Gj. The correlator (£ m ,£„) = 7<5m,n (where for the 
cosine phase factors 7 = 1/2), and we obtain 

<SGi,Ea ti >=< (Gi, &)(Ei, &) >=7 E i G i &,^) 

(34) 

The scalar product of the normalized positively-defined 
vectors in the brackets is of the order of unity, and one 
comes to the estimate for the average, Eq.||, < ... >= 
CSGE e i (according to definitions given in front of Eq.7, 
E = E el and G = SG. 

Representing vector G as a sum of components "par- 
allel" and "normal" to the vector E: G = + G± one 
has (G B ,E) = (G,E), (Gj_,E) = 0, and, finally, 

G, = a^E = cf E 

E 2 E 

Decomposing in the same way the random vector 
on the components "parallel" and "normal" to E (fg 
and and taking into account that 

<SG> E =< (G B + G_ L ,£ S + £ L ) > E , 



and < (Gj_,£j_) >= we finally have 

< SG > E =< (G e ,Ze) >= <^(E,f) = C 6 S E 

hi hi 

which corresponds to Eq.^. It means that SG has a linear 
regression with respect to E e i . 
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